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ABSTRACT 

We report observations of the remnant of Supernova 1987A with the High Resolution Camera 
(HRC) onboard the Chandra X-ray Observatory. A direct image from the HRC resolves the annular 
structure of the X-ray remnant, confirming the morphology previously inferred by deconvolution of 
lower resolution data from the Advanced CCD Imaging Spectrometer. Detailed spatial modeling shows 
that the a thin ring plus a thin shell gives statistically the best description of the overall remnant 
structure, and suggests an outer radius 0.96"±0.05"±0.03" for the X-ray-emitting region, with the two 
uncertainties corresponding to the statistical and systematic errors, respectively. This is very similar 
to the radius determined by a similar modeling technique for the radio shell at a comparable epoch, 
in contrast to previous claims that the remnant is 10-15% smaller at X-rays than in the radio band. 
The HRC observations put a flux limit of 0.010 ctss"^ (99% confidence level, 0.08-10 keV range) on 
any compact source at the remnant center. Assuming the same foreground neutral hydrogen column 
density as towards the remnant, this allows us to rule out an unobscured neutron star with surface 
temperature T°° > 2.5 MK observed at infinity, a bright pulsar wind nebula or a magnetar. 

Subject headings: circumstellar matter — shock waves — supernovae: individual (SN 1987A) — 
supernova remnants — X-rays: general — stars: neutron 



1. INTRODUCTION 

The core-collapse supernova (SN) 1987A in the Large 
Magellanic Cloud was the brightest SN observed since 
the invention of modern telescopes, providing the best 
opportunity to study the last evolutionary stage of a mas- 
sive star. Optical observations have revealed a triple-ring 
nebula centered on the explosion, consisting of an inner 
equato rial ring of radius 0.8 1 "-0.86"and t wo la rger outer 
rings (jBurrows et all Il995t iPlait et al.l I1995D . These 
are part of the hourglass-shaped circumstellar medium 
(GSM) ejected by the progen itor 20,000 years ago (see 
iMorris fc Podsiadlowskill2007f l. Since early 2004, the SN 
blast wav e has begun to encount er the main body of in- 
ner ring (jPark et al.ll2005L l2006h . This 'big crash' has 
led to a drastic increase in soft X-ray emission that orig- 
inates from the optically thin thermal plasma behind 
the shock. The X-ray spectrum is well-described by a 
two-component plane-parallel shock model in nonequi- 
librium ionization, with plasma temperatures kT ~ 2 
and 0.3 keV, correspo nding to the fast and d ecelerated 
shocks, respectively (|Park et al.l 12004 120061 ). As the 
blast wave propagates into the dense GSM, the soft X- 
ray emission traces the evolution of the forward shock, 
probing the GSM structure and its density profile. 

The nearness of SN 1987A (51.4 kpc) allows us to re- 
solve a supernova remnant (SNR) morphology at a very 
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young age. In X-rays, this is only possible with the 
Chandra X-ray Observatory, the highest resolution X-ray 
telescope compared to any previous, current and even 
planned future X-ray missions. Previous Chandra ob- 
servations of SNR 1987A were all carried out by the 
Advanced CCD Imaging Spectrometer (ACIS). However, 
the ACIS detector has a pixel size 0.492", comparable to 
the full width at half maximum (FWHM) of the mirror's 
on-axis point-spread function (PSF). Although the effec- 
tive resolution can be slightly improved by the dither- 
ing of the spac ecraft and by applyin g a sub-pixel imag- 
ing algorithm (|Tsunemi et al.ll2001| ). image deconvolu- 
tion is sti ll required to fully re solve the remnant struc- 
ture (e.g. iBurrows et al] 120001 ). To determine whether 
artifacts are introduced by the complicated non-linear 
reconstruction process, the ACIS results need to be com- 
pared with higher resolution direct X-ray images. An 
analogy can be drawn from the radio imaging campaign 
of SNR 1987A. Since 1992, Austraha Telescope Compact 
Array (ATCA) observations at 9 GHz have revealed de- 
tailed structure of the radio shell using the super-resolved 
technique, but it was not until the upgrade of the ATCA 
in 2003 that the higher resolution diffraction-limited im- 
ages at 18 GHz provide d a direct i mage of the radio mor- 
phology (see review bv lGaensler et al . 20071). 

The High Resolution Camera (HRC) onboard Chan- 
dra, consisting of two microchannel plate detectors, offers 
smaller electronic readout pixels (0.13175") than ACIS. 
These better sample the PSF, providing a more straight- 
forward imaging process without the need for deconvolu- 
tion. Despite a smaller effective area for HRC than ACIS 
and the lack of any spectral resolution, the SNR is now 
very bright in X-rays below 2 keV, at which the HRC has 
good sensitivity, making it an ideal instrument for mor- 
phological studies. In this Letter, we report a detailed 
analysis of the first HRC observation of SNR 1987A. 
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Fig. 1. — (a) Raw HRC data of SNR 1987A taken on 2008 Apr 28-29. The image was binned into the HRC detector pixel and centered 
at RA=05'^35""28=, Dec=-69°16'11.2" {J2000). (b) Raw ACIS data taken on 2008 Jan 9. The image is in 0.3-8keV energy range and 
was binned into the ACIS detector pixel, (c) Super-resolved radio image at 9 GHz taken by ATCA on 2008 Apr 23 {Ng et al. 2008). (d) 
Deconvolved HRC imago using the dataset shown in panel (a), (e) Deconvolved ACIS image using the dataset shown in panel (b), in 
0.3-8 keV energy range (Racusin et al. 2009). (f) HRC data in panel (a) smoothed to 0.4" to match the resolution of the radio image in 
panel (c). All panels are on the same spatial scale. 



2. OBSERVATION AND DATA REDUCTION 

Our Chandra observations were carried out on 2008 
Apr 28-29 (day 7736 since the SN explosion) using the 
HRC-I detector with a total exposure of 46 ks (ObsIDs 
9085 and 9851). All the data reduction was carried out 
using CIAO 4.1^ with CALDB 4.1.1. In addition to stan- 
dard data pro cessing, we applied a screening algorithm 
(jMurrav et al.l 12000) to reject the 'bad events' including 
non- X-ray background and mislocated events, thus im- 
proving the data quality and boosting the signal-to-noise 
ratio. There were no strong background flares during the 
exposure; hence all data were included in the analysis. 

3. SPATIAL ANALYSIS AND RESULTS 

The HRC image of SNR 1987A is shown in Figure [TJi. 
This reveals a ring-like morphology for the remnant with 
a subluminous region (a 'hole') at the center. The X-ray 
emission peaks in the northeast, and the north rim is gen- 
erally brighter than the south. A comparison to the raw 
ACIS image at a similar epoch (2008 Jan 9) is shown in 
Figure [TJd. The HRC image resolves, for the first time, 
the remnant structure with direct X-ray imaging, and 
reveals the central hole unambiguously. The entire rem- 
nant has (3.0 ± 0.02) x lO'' cts in the full energy range 
(0.08-10 keV) of the HRC and the background is neghgi- 

® jhttp : //cxc ■ harvard ■ edu/ciao7] 



ble (~ 0.1%), implying a count rate of 0.66±0.01ctss^^. 

To further improve the resolution of the HRC image, 
we applied a s imilar deconvolution process as used by 
IBurrows et all ( 2000). The data were binned to a quarter 
of the detector pixel in each dimension, then deconvolved 
using the Lucy-Richardson algorithm (|Lucv.,1974) with 
a ChaRT/MARX'^-simulated PSF at IkeV, where both 
the source spectrum and the detector response peak.^ 
The resulting image after 30 iterations was smoothed to 
0.2" resolution and is shown in Figure [TJi. As a com- 
parison, we also showed in Figu re [it the deconvolution 
ACIS image from lRacusin et all (2009). In both images, 
the remnant exhibits a similar annular structure with 
sharp edges, and has semi-major and semi-minor axes 
of 0.8" and 0.6", respectively. There are some very mi- 
nor difference in the extent of the bright regions between 
the two images, but the brightness variation is mostly 
continuous and smooth around the ring, and there is no 
obvious point source at the remnant center. 

Figure [TJ: shows the super-resolved radio image of 
SNR 1987A at 9 GH z, taken by the ATCA at a similar 
epoch (2008 Apr 23: INg et al.1l2008D . To obtain a direct 
comparison between the radio and X-ray images, we have 
regridded and smoothed the HRC data to the same reso- 
lution of 0.4" FWHM, and the resulting image is shown 



^ http://cxc.harvard.edu/chart/ 
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Fig. 2. — (a) Surface brightness profiles of SNR 1987A extracted from elliptical annulus regions (see text), (b) Slices through the remnant 
at different position angles (north through east). The offset is radial from the center and is positive toward north and/or west. In both 
plots, the raw and deconvolved HRC images are shown by the histograms and hnes, respectively. 



in Figure [Tf. The radio shell consists of two main lobes: 
a bright one in the east and a fainter one in southwest; 
it has a similar size as the X-ray ring and both emissions 
peak in the east. Also, there is some hint of an X-ray 
counterpart for the southwestern radio lobe. However, 
except these similarities, the radio and X-ray emissions 
exhibit substantially different morphologies. While the 
former is mainly contributed by the eastern lobe, the lat- 
ter is brighest along the north rim, and neither the overall 
bilateral pattern nor the strong east-west asymmetry in 
the radio image are observed in the HRC data. 

We extracted surface brightness profiles for both raw 
and deconvolved HRC images using elliptical annulus re- 
gions with aspect ratios fi xed by the system's orientation 
(43.4° to the line of sight: [Pimlf200l . and the results are 
plotted in Figure [5^. The X-ray emission is very faint 
at the remnant center and peaks at about 0.8" radius 
with relatively sharp boundaries. The profile of the de- 
convolved image indicates a width of FWHM 0.5" for the 
ring, but we note that the actual width is likely smaller, 
since the deconvolved image has been smoothed. Fig- 
ure [^b shows slices through the remnant along different 
position angles, illustrating a high degree of symmetry 
along the east-west direction. Wc found that counts in 
the eastern and western halves differ only by 5% ± 3%. 
In contrast, the north-south asymmetry is more signifi- 
cant (25% dz 7%). It is partly due to the Hght travel time 
effect: the rapid brightening of the remnant makes the 
north rim, which is currently 1.5 light year closer to the 
Earth than the south rim, appear brighter. Based on 
the ACIS lightcurve, we estimated that this effect could 
result in at most 15% count difference between the north- 
ern and southern halves. However, it is not enough to 
account for the observed value; therefore, we conclude 
that the rest is intrinsic. 

Motivated by the images, we have carried out detailed 
spatial modeling to capture the characteristic scales of 
the SNR. Simple models were generated in 3D, corrected 
for light travel time effects and projected onto the image 
plane according to the viewing geometry. The resulting 
images were then convolved with the PSF and fitted to 
the raw HRC data show n in Figurc^i using the pro- 
cedure described by iNg fc Romanil (120081 ) with iGehreld 



(jHH) X statistic. Since we a re mostly intere sted in 
the global geometry, we followed iNgeTall (pOOl to ac- 
count for the overall asymmetry using a linear brightness 
gradient with arbitrary position angle in the equatorial 
plane. It is worth noting th at as compared to the recent 
ACIS study (jl acusm "etall [20091. our fitting procedure 
involves no resampling or deconvolution processes that 
may degrade the data or introduce artifacts in the im- 
ages. Therefore, it can provide robust measurements. 

We have tried fitting various mo dels, including a simple 
ring, an equatorial belt torus (see INg et al.ll2008l ). and a 
ring plus a shell. The best-fit parameters are listed in Ta- 
ble [T] with the corresponding images shown in Figure [H 
Note that the formal reduced values for the fits are 
likely underestimated, since the number of independent 
pixels in the image we are fitting is much larger than the 
source area. Better estimates of the reduced values 
could be made using the source area of ^ 500 pixels. As 
shown in the table, this gives more reasonable results. 
The uncertainties quoted in Table [T] are statistical errors 
with 90% co nfidence level, estirn ated from Monte Carlo 
simulations (|Ng fc Romani|[2b08[ ) . We also followed their 
procedure to quantify systematic errors due to features 
not captured by the simple models, such as small-scale 
brightness fluctuations around the ring. This suggests 
a systematic uncertainty of ~ 3% for the radius. As 
indicated in Table [H different fits give similar radii of 
~ 0.8", consistent with the plot in Fig. [5^. We found 
that a thin ring fits poorly to the data, but varying the 
ring's thickness significantly improves the result. Our 
thick ring model is essentially s ame as the sim ple torus 
used in the ACIS study of Racusin et al.l ()2009f ). and the 
radius and thickness we obtained are in good agreement 
with their results. Our equatorial belt torus model of- 
fers a better fit than a ring, and suggests a half-opening 
angle ~ 26 °, slightl y smaller than that of the radio shell 
(|Ng et al.i r2008: P otter et al.ll2009( ). However, the char- 
acteristic two-lobe structure of the equatorial torus does 
not match the data (Fig. [Hb)- In terms of the statis- 
tics, our best model is given by a thin shell plus a thin 
ring. As shown in Figs. & Os, this model captures 
the overall remnant morphology well and successfully re- 
produces the observed ellipticity. The fit suggests 18% 
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TABLE 1 

Best-fit Spatial Models to 2008 Apr, 28-29 Chandra HRC Image of SNR 1987A 



Model Ring radius (") Shell radius (") Thickness {") dof xi 

Thin ring 0.83 ± 0.01 • ■ ■ 0.04=^ 1218 4090 2.4 

Thick ring 0.75 ±0.02 ■ ■■ 0.56 ± 0.08 942 4089 1.9 

Torus'^ ■■■ 0.82 ±0.02 0.39 ± 0.02 884 4088 1.8 

Ring+shelP 0.82 ±0.03 0.94 ± 0.05 0.05^^ 709 4088 1.4 



Note. — The uncertainties quoted are statistical errors at 90% confidence level. 
^ The formal number of degrees of freedom (dof) depends on tlie image size we 
are fitting, and they are likely overestimated (see text). 

^ The reduced values, Xvj ^Lve estimated with dof=500, as we argued in the 
text. 

Fixed at 5% of the corresponding radii during the fit. 
The best- fit torus has a half-opening angle 26° ± 3° . 
° The fit suggest 18% more counts in the ring than the shell component. 



more counts in the ring component than the shell. While 
it is tempting to associate the latter with faster shocks 
at a higher latitude, the physical scenario is likely to be 
more complicated and the parameters in the fits could 
be highly degenerate. Further observations are needed 
to confirm if a two-componen t spatial model is required. 

Similar to INg et all ([2008), we added a point source 
at the center of the best-fit model to estimate the flux 
limit from any possible central object. We adjusted 
the point source flux and re-fit the model until the 
value exceeded a certain level. This gives a de- 
tection limit O.OfOctss"^ at 99% confidence level over 
the energy range 0.08-10 keV. For a = f3km ra- 
dius^ neutron star 51.4kpc away with a foreground neu- 
tral hydrogen colu mn density N-^ — 1.3 x 10^^ cm^^ 
(jZhekov et al.ll2009t ) , the flux limit corresponds to a sur- 
face temperature T°° = 2.50 MK, or a bolometric lumi- 
nosity Lj^j = 4.7 X 10^* ergs s^^. More physical estimates 
using neutron star atmosphere models (e.g. Zavli n et al.l 
Il996| i suggest similar values. If the radiation from the 
central source is nonthermal with a typical powerlaw in- 
dex r = 1.5, then the limit converts to a luminosity of 
7.0 X 10^'*ergss~^ in 2-lOkeV energy range. To have 
a fair comparison with previously published values, we 
followed Park et al. (2004) to ignore the statistical fluc- 
tuations of the underlying remnant, and repeated the 
above exercise. This gives a 90% confidence level de- 
tection limit 3.3 x lO^^ergss^^ in the 2-lOkeV range, 
slightly higher tha n the ACIS hmit of 1.5 x lO'^^ergss"^ 
(jPark et al.l 120041 ) . likely due to the brightening of the 
remnant, but more stringent than the XMM-Newton 
limit 5 X 10'^^ ergs s^^ o btained from spectral analysis 
(|Shtvkovskiv et al.l[2005h . 

4. DISCUSSION 

4.1. X-ray Morphology of SNR 1987A 

The HRC observations of SNR 1987A have provided a 
direct X-ray image that confirms the remnant structure 
as suggested by previou s ACIS obse rvations using the 
deconvolution technique (jRacusin et al. 2009, and refer- 
ences therein). As the image shows, the entire X-ray 
ring has been lit up, and its size and ellipticity are in 
good accord with the optical inner ring. This gives di- 
rect evidence that the blast wave has already encountered 

^ This is the redshifted value as observed at infinity, same for 
the temperature and luminosity below. 



the dense CSM all around the ring and is now interact- 
ing with the main body of the ring. The deconvolved 
HRC and ACIS images are nearly identical, except for 
some very minor differences in the brightness distribu- 
tion, which could be attributed t o the dif f erent energy re- 
sponse of the instruments. Zhck ov et al.l (j2009l ) reported 
a strong energy dependence of the east-west asymmetry 
in the ACIS data. The eastern half is dominated by hot- 
ter plasmas behind fast shocks, therefore it is brighter at 
higher energies. Conversely, the western half has cooler 
plasmas, slower shocks and is brighter at lower energies. 
These two contributions could possibly cancel out in the 
HRC image, resulting in the high degree of symmetry 
observed along the east-west direction. We are unable 
to verify this because the HRC data do not possess any 
spectral information. As the shocks further propagate 
and decelerate in the dense CSM, we expect the degree 
of symmetry to vary in future observations. 

The poor spatial correlation between the soft X-ray 
and radio morphologies suggests that they may origi- 
nate from different physical processes. While the X- 
rays trace the SN blast wave, the radio emission is 
believe d to originate between t he forward and reverse 
shocks (jManche stcr et al.l 120051 : INg et al.l l2008 f) . Given 
the young age of S NR 1987A, both shocks should have 
similar radii (e.g. iTruelove fc McKed 1 1999( 1. hence, one 
would expect the X-ray and radio shells to have a com- 
parable size. A direct comparison between the HRC 
and radio images in Figure [1] seems to support this pic- 
ture. Nonetheless, previously reported radii for the X- 
ray shell h ave always been 10-15% smaller th an for the 
radio shell (|Ng et al.ll200llRacusin et al.ll2009f ). The dis- 
crepancy could Ije due to different rn easuring techniques 
(jCaensler et al.|[2007t iRacusin et al.|[2009i) . or due to sys- 
tematic errors in the measurements. A better way to 
compare the sizes is from the outer edge of the shells. 
This can avoid complication by the inner ring or projec- 
tion effects. After accounting for the shell's thickness, 
our best-fit model of the HRC image has an outer radius 
of 0.96". Similar values are obtained from other mod- 
els (see Tab le [It , and it also seems consistent with the 
ACIS study (jRacusin et al. 2009). All these indicate that 
our estimate is robust, thus, we conclude that the X-ray 
shell has an outer radius of 0.96" ± 0.05" ± 0.03", with 
the uncertainties corresponding to statistical (90% con- 
fidence) and systematic errors, respectively. This is very 
similar to the outer radius 0.92" ± 0.06" of the radio 
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Fig. 3. — Best-fit spatial models for the HRC image of SNR 1987A: (a) a thick ring, (b) a torus, (c) a thin ring plus a thin shell. All 
images are convolved with the instrument PSF. (d) The HRC data, same as Fig. [3ji. (e) MARX simulation of the best-fit ring plus shell 
model from (c). (f) Difference image between panels (d) and (c), showing the residuals of the fit. All panels except (f) are on the same 
color range. 



shell (jPotter et al.ll2009l ). thus confirming the physical 
picture discussed above. Although these values appear 
to be slight ly larger tha n the opt ical inner ring of radius 
0.81"-0.86" ((Burrows et al. 1995: IPlait eFall flQQS). they 
unlikely represent the radius of the transmitted shock in- 
side the dense CSM, since which should have a very slow 
propagation s peed according to some hydrodynamic sim- 
ulations (e.g. lDwarkadasll2007f) . Therefore, we believe 
that our results could correspond to the projected size of 
fast shocks traveling in a lower density environment at 
higher latitudes. 

4.2. X-ray Limit on a Central Source 

In the absence of fast cooling mechanisms, stan- 
dard theories suggest a bolometric luminosity > 
5.0 X 10'^'* ergs s~^ (corresponding to a surface temper- 
ature T°° ^ 2.54 M K) for a 21-year-old neutron sta r 
(lYakovlev fc Pethickl (2004t IShternin fc YakovlevI [200l . 
Our flux limit rejects any unobscur ed neutron stars hot - 
ter than this at 99% confidence. [Graves et al.l ()2005D 
argued that a neutron star obscured by dust is unlikely, 
since the IR flux is consistent with the radioactive decay 
of "^^Ti in the ejecta and no other energy source is re- 
quired. The s ame conclusion also holds for the recent IR 
observations (jBouchet et al.l [20061 ). On the other hand, 
for a neutron star with short thermal relaxation time due 
to superfluidity, quark matter, or anomalously high ther- 
mal conductivity of the crust, the thermal emission from 
the surface could have already been too low to detect 



at an age of 21 (|Shternin fc YakovlevI 120081 : iChan et all 
l2009l) . 

X-ray emission from a young pulsar is often dominated 
by nonthermal radiation from the magnetosphere and 
from the surrounding pulsar wind nebula (PWN). The 
nonthermal flux limit we derived rules out energetic pul- 
sars such as the Crab or PSR B0540-69, but a fainter 
pulsar/PWN s ystem similar to 3C 58 wou ld still escape 
detection (see iKargaltsev fc Pavlov! I2008f ). These au- 
thors also gave a general correlation between nonther- 
mal luminosity and spin-down power E of young pulsars, 
which places a conservative limit _E < 3 x 10'^'^ ergs s~^ 
for the unseen pulsar in SNR 1987A. Based on the 
non-detection, lOgelman fc Alpail (|2004l ) argued that the 
pulsar, if it exists, should have either a very strong 
or very weak magnetic field. While the former case 
would put the star into the magnetar regime, our de- 
tection limit is incompatible with the X-ray luminosi- 
ties of nearly all known rn a gnetars, even in quiescence 
(jWoods fc ThompsonI [200l iMereghettil 12008^. For the 
weakly magnetized neutron star case, one candidate is 
a central compact object (CCO), which is a soft X-ray 
source found inside a SNR, and could possibly be a neu- 
tr on star born with a weak magnetic field (see review 
bv lGotthelf fc Halpern|[2008f ). It has been proposed that 
SNR 1987A may ha rbor a CCO (e.g. Manchester 200'^; 
IHalpern et al.l [20071 ): however, the current observations 
do not provide any useful constraints on this scenario, 
unless CCOs have similar cooling properties as discussed 
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above. Finally, we note that [Graves et alj (|2005[ ) con- 
sidered different accretion models and used the UV flux 
limit to rule out spherical accretion and a slim disk, 
where the disk thickness is non-negligible. They also 
placed stringent constraints on thin-disk models. Al- 
though it is beyond the scope of this paper to repro- 
duce their modeling, given the similar order of our X-ray 
limit, we believe that our results unlikely improve their 
constraints substantially. 

5. CONCLUSION 

We have presented a detailed analysis of Chandra HRC 
observations of supernova remnant 1987A. The remnant 
has a similar size as the optical inner ring, and its mor- 
phology is statistically best-fitted by a thin ring plus a 
thin shell. This provides direct evidence that the dense 



inner ring of circumstellar medium has been overtaken 
the supernova blast wave. Our spatial modeling also 
indicates a similar size between the X-ray- and radio- 
emitting regions, confirming the picture that the super- 
nova forward and reverse shocks are closely located. The 
X-ray flux limit on any possible central source rejects a 
young neutron star, unless it has some fast cooling mech- 
anisms or is obscured by a small accretion disk. 
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